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We examine the state of statistical equilibrium attained by a uniformly forced condensable sub- 
stance subjected to advection in a periodic domain. In particular, we examine the probability density 
function (PDF) of the condensable substance in the limit of rapid condensation. The constraints 
imposed by this limit are pointed out and are shown to result in a PDF — whenever the advect- 
ing velocity field admits a diffusive representation — that features a peak at small values, decays 
exponentially and terminates in a rapid "roll-off" near saturation. Possible physical implications of 
this feature as compared to a PDF which continues to decay slowly are pointed out. A set of sim- 
ple numerical exercises which employ lattice maps for purposes of advection are performed to test 
these features. Despite the simplicity of the model, the derived PDF is seen to compare favourably 
with PDF's constructed from isentropic specific humidity data. Further, structure functions asso- 
ciated with the condensable field are seen to scale anomalously with near saturation of the scaling 
exponents for high moments — a feature which agrees with studies of high resolution aircraft data. 

PACS numbers: PACS number 47.52.+j, 05.45.-a 



I. INTRODUCTION 

Advection- Diffusion- Condensation (ADC) is a variant of the familiar passive advection-diffusion problem Q], wherein 
apart from molecular diffusion a tracer is subject to an additional sink associated with the process of condensation. A 
study of the interplay of theseprocesses is motivated by the need to understand the large scale distribution of water 
vapor in the troposphere Indeed, the prominent role played by water vapor in various problems related 

the Earth's climate [j| — for example, via the water vapor feedback in the greenhouse effect |(J — makes this an 
important issue. In this context, ADC serves as an idealized problem for water vapor in the troposphere 0. The 
mixing ratio q(x, t) of a scalar tracer subject to advection, diffusion and condensation is governed by 



■^ + (u-V)q = nX7 2 q + S(q,q s ) + F (1) 

here u and F are the advecting velocity field and the forcing respectively and rj is the diffusivity of the condensable 
substance. In the present work we shall assume the domain to be spatially bi-periodic and restrict our attention to 
the nondiffusivc limit, r\ = 0. The problem thus consists of the advection of a passive tracer supplied by the source 
F, and removed by a sink associated with condensation. This sink is represented as 

S = --[q - q s (x}] if q > q s 

T 

= if q<q s (2) 

where q s (x) is the saturation mixing ratio — a prescribed function. Hence with r\ = 0, the ADC problem is closely 
related to the interaction of smooth advection with linear damping |8|. The ADC problem also bears a similarity 
recent studies of "active" processes when coupled with fluid advection, such as the evolution of chemically [9j or 
biologically ^3] active substances and the issue of phase separation in immiscible fluids . 

Our aim is to examine the state of statistical equilibrium attained by and |J2J in the limit r — > or that 
of rapid condensation. This limit implies that upon advection if q(x,t) > q s {x), then the parcel's mixing ratio is 
instantaneously reset to the saturation mixing ratio at that location. Of course, two immediately apparent opposing 
long time limits of Q and 101 ar e : (i) Uniform forcing with u = q(x, t) — > q s (x) and (ii) A sufficiently mixing 
flow with F = => q(x.t) — > mhLyg[q s (x)]. Hence, to achieve a non-trivial state of statistical equilibrium we require 
both u, F j= 0. 
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In the absence of diabatic effects, parcel motion in the midlatitude troposphere is restricted to two-dimensional sur- 
faces of constant potential temperature, i.e. to isentropic surfaces Our objective is to understand the probability 
density function (PDF) of the water vapor mixing ratio along these midlatitude isentropic surfaces O n average 
the thermal structure of the midlatitude troposphere is such that as one progresses polewards, the temperature along 
isentropic surfaces decreases in a fairly linear manner. Given the sensitivity of the saturation mixing ratio (via the 
Clausius-Clapeyron relation) to the temperature, following Q we take q s to vary exponentially with y. As we are in 
a periodic domain, q s (x) = q s (y) = exp(— a\y — ^|) (0 < y < L) where a > and increasing a yields progressively 
steeper profiles which fall off symetrically from y = This yields an idealized model problem which retains the 
essence of the much more complex atmospheric problem that motivates our work. Even though the real atmosphere 
does not conform exactly to the idealizations, progress can be made through a detailed solution of this model problem. 



II. THE PDF IN THE LIMIT OF RAPID CONDENSATION 



A. General formulation and boundary conditions 



Before we proceed to the PDF in the general case, consider the limit alluded to in the Introduction. Specifically, 
u = with constant forcing. On one hand this case is quite straightforward but on the other it exhibits some 
pathologies inherent in dealing with the PDF. From QJ at a location xq, yo f° r t> t\ = q s (xo, Uo)/F we have 



q(x ,y ,t) = q s (x ,y ) + Ft[1 - cxp(— — )] (3) 

T 

As © shows, in general the domain will become supersaturated. But for rapid condensation, as claimed, 
l{xo, 2/0j t) Qs(xo, yo) V(xo, yo)- Note that the limit of saturation is approached from above. Physically this implies 
that as t — ► 0, we simultaneously have q(x,y,t) — > q s {x,y) such that the condensation sink balances the forcing. 

To estimate the PDF (denoted by P = P(x,y,q,t)), we decompose the stationary solution as 
P(x,y,q) = P\(x,y)Fi(q\x,y) where F\(q\x,y) is the conditional PDF of q conditioned on (x,y). From the 
preceeding discussion, F\{q\x,y) = S[q — g s (a;,y)]. As u = 0, for an unbiased estimate it is natural to take P\(x,y) 
to be uniform resulting in P(x,y,q) = S[q — q s (x,y)]. In essence we end up with a S function supported on the 
fixed point of the dynamical system Q with u = 14] . In this particular case, the system does not have any 
other invariant PDF and the S function PDF is meaningful (as will be seen later), but in general to avoid such 
complications we restrict our attention to functions that are reasonably smooth. 

Proceeding to the general case pq. the Liouville or transport equation satisfied by the PDF is 



dP(x,y,q,t) | d(u t P) | d[(F + S)P] _ Q 
dt dxi dq 

with / ( P(x,y,q,t) dx dq = IVt (4) 

Jx JO 

as in general, i.e. outside the limit of rapid condensation, < q(x,y,t) < oo and apart from periodicity in x we 
require P(x, y, oo, t) = 0. 

In (x,y,q) space the above can be viewed as resulting from a probability current UP where U = [u,v, (F + S)]. 
An immediate consequence of rapid condensation is a restriction on q i.e. q(x,y,t) < q a (x,y). Till now our choice of 
forcing has been completely arbitrary, at this stage we restrict ourselves to a forcing which is positive definite and also 
satisfies ^ = 0. This implies, whatever the initial condition on q is, after a finite time mm.{q) = mm(q s ). Hence, 
min(g s ) < q(x,y,t) < q s (x,y) and using (0), @J reduces to 
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dP(x,y,q,t) | djujP) d(FP) =Q 
dt dxi dq 

r r ls(x,y) 

with / / P(x,y,q,t) dx dq = 1 Vt (5) 

J x J min (g s ) 

and, as noted, apart from the restriction on the forcing function, we require P to be relatively smooth. Also, note 
that now the domain is periodic in (x,y) and is bounded by the surfaces q — min(g s ) and q — q s (x,y). Further, 
periodicity and the imposition of rapid condensation imply [J7P]|aD = and in essence (J3J) represents the evolution 
of a PDF via an incompressible flow within a bounded domain with the normalization constraint implicity providing 
the required boundary condition. Interestingly, similar integral formulations for the PDF arise in studies of spiking 
neurons (see for example Fusi & Mattia [l5| and the references therein). 



B. Solutions for present saturation profile 



In the present case, as the saturation mixing ratio is purely a function of y, if the forcing is also taken to be of 
the form F = F(y) (in fact we will focus on the uniformly forced case) — it is reasonable to look for solutions 
which are independent of x. Of course, it is the realization of statistical equilibrium without the presence of gradient 
fields (as ij — 0) that makes progress possible in the present case [2]). This circumvents the need to estimate condi- 
tional expectations of the scalar dissipation (or diffusion) which complicate advection-diffusion problems ^1,^3' E3- 



1. The uniform case 



Inspecting (J5J it is evident that without any further assumptions, as per the usual situation involving the evolution 
of a PDF via an incompressible flow, a valid stationary solution is that of a uniform distribution in (y,q). To obtain 
the PDF of q from this uniform distribution we estimate the probability that q < Q as a function of Q (i.e. Q 
represents the sample space variable corresponding to q). Utilizing the rapid condensation normalization constraint, 
P(q,y) has to be integrated with respect to the saturation curve q s (y). Therefore 

Pr{q <Q) = (6) 
it 

Here ii, I2 correspond to the regions y <Y and y > Y where Y = Z(Q) (Z being the inverse of q s (y)). In the present 
case, as q s is symmetric about y = ^ we need only consider half of the domain (L/2 < y < L) with q s — exp(— ay), 
i.e. Z(Q) = — \og(Q)/a. Specifically, 

z(Q) r Q 



h= / P(q,y) dqdy 

JL/2 JO-? 



'L/2 JQ-2 
r L r <ls(y) 

h= / P(q,y) dqdy (7) 

J Z(Q) JQ 2 

where Q2 = min(g s ). I t in © is the same as I2 but with L/2 as the lower limit of integration in the outer integral, 
which ofcourse is nothing but the normalization constraint in JHJl- Hence, 

PDF(0) . £l±M (8) 

Substituting P(y,q) = const, in Q and © yields PDF(Q) ~ Z(Q) i.e. PDF(Q) ~ log(Q^) which can be seen in 
the upper panel of Fig. In essence, for a uniform distribution the constraint of rapid condensation forces the 

PDF to reflect the saturation profile. 
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2. The eddy- diffusion case 



We now assume that the the effect of the fluctuating velocity in the Liouville equation admits a diffusive 
representation — say by an eddy diffusivity K e . This would be the case if the parcel trajectories consisted of 
independent Brownian motion, for example. Admittedly, this is a fairly severe assumption in that mixing by multiple 
scale velocity fields rarely follows a simple diffusive prescription |l9j . But from a tropospheric viewpoint, it is known 
that the meridional (y - direction) Lagrangian eddy-velocity correlation function decays quite rapidly (on the order 
of a few days) |2(lj — hence, in this context the assumption of an eddy diffusivity may not completely unreasonable. 
Note that the eddy diffusivity we refer to here is an eddy diffusivity applied to probability evolution in (y, q, f) space. 
This is not the same as characterizing mixing by an eddy diffusivity to an evolution equation for a coarse-grained q 
in (y, t) space. Indeed, in Q it is shown that the Brownian model yields different coarse-grained q statistics than the 
mean-field diffusivity model; it is suggested further that the Brownian motion model constitutes a minimal model for 
investigation of the interplay of transport processes with a nonlinear sink term such as condensation. In that sense, 
the Brownian model is worth of study in and of itself. 



With the eddy-diffusivity JSJ yields 



dP d 2 P „, ,dP 



K e^+F(y)— =0 (9) 



For uniform forcing F(y) = S (a constant), as P > and is periodic in y, by inspection the ansatz 



1 K 

P(q,y) ~ exp(-fc 2 <7)[cos(Ay) + sin(Ay)] ; where A = -,k 2 = \ 2 -f (10) 

4 o 

satisfies a stationary form of JHJ). Indeed, for a quantitative estimate we would need to fix the constants that arise 
in HI U| l via the normalization constraint. But for a quantitative estimate substituting P(q,y) in Q and setting 
dG/dy = [cos(Ay) + sin(Ay)] for notational simplicity, yields 



PDF(Q) ~ cxp(-fc 2 Q) [G(Z(Q)) - G(0)] (11) 

Substituting for G(y) 

PDF(Q) ~ 6Xp( ~ fc2Q) {sin[AZ(Q)] - cos[AZ(Q)] + 1} (12) 
A 

and finally substituting for Z(Q) 

PDF(Q) ~ CXp( A fc2Q) {sin[log(Q^ )] - cos[log(Q^ )] + 1} (13) 

The dependence of the PDF on is shown in the lower panel of Fig. JIJ. In all cases, the PDF has a peak for small 
q decreases exponentially for intermediate values of q (the slope increases with i | L ) and then rolls-off as q — > max (q s ). 

Comparing this with the uniform case (upper panel of Fig. (JTJ), it is evident that the form of the PDF is still 
controlled by the saturation mixing ratio profile. But the details, such as the slope of the PDF for intermediate 
values of Q, are dependent on the eddy-diffusivity and forcing. Also, it is worth noting that the PDF is in 
marked contrast to what one would obtain for a fully saturated domain. For example in the saturated case we had 
P{q.iV) = $(q ~ Qs(y)), this yields PDF(Q sat ) ~ —dF(Q)/dQ which implies a power law (Q _1 , with no "roll-off" for 
high Q) when q s — exp(— ay). At first sight the roll-off in l|13f) appears to be an obvious result of enforcing rapid 
condensation, but note that the saturated case considered above also conforms with the rapid condensation limit yet it 
yields a very different PDF. This is illustrated more clearly by choosing q s (y) to be a linearly decreasing function of y, 
now the saturated PDF is uniform, whereas i|12|) again yields a slowly decaying PDF with a smooth roll-off at large q. 
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With regards to water vapor, given the logarithmic dependence of the infrared cooling on the water vapor mixing 
ratio, it follows that fluctuations in the mixing ratio actually increase infrared cooling |7j. Hence a slowly decaying 
PDF, by increasing the probablity of encountering a large fluctuation (as compared to a normal PDF), increases 
the efficiency of radiative cooling to space. In this regard, the rapid roll-off of the tail of the PDF could play 
a strong role in that it sharply decreases the possibility of very large fluctuations — for example, in comparison 
to a advection-diffusion model of water vapor distribution where one would encounter exponential tails under 
homogeneous forcing |2l| — and would imply a reduction in the infrared cooling or in other words an increase in the 
temperature of the surface to maintain energy balance. Quantifying this effect by using a full-fledged atmospheric 
radiation code is a project we hope to pursue in the near future. 



III. NUMERICAL SIMULATIONS 

A primary assumption in our derivation was the use of an "eddy diffusivity" to represent the effect of an advecting 
velocity field. As mentioned, in most flows of interest the velocity field is expected to be quite coherent and might 
not conform to a diffusive representation. To test the stringency of this assumption we numerically simulate the ADC 
system by employing a lattice map for purposes of smooth large scale advective mixing (2^ , |2^| . The velocity field is 

u(x, y, t) = f(t) Ai s\n{Biy + p n ) 
v(x, y, t) = (l- f(t)) A 2 sm(B 2 x + q n ) 

(14) 

where A 1 = 0.75, A 2 = 0.5, B 1 = B 2 = 2.5, /(*) is I for nT < t < (n + l)T/2 and for (n + l)T/2 < t < (n + 1)T. 
PmQn (€ [0, 27r]) are random numbers selected at the beginning of each iteration, i.e., for each period T. Advection 
is implemented via sequential integer shifts in the x and y directions on a square lattice (see |22( for details). Apart 
from its numerical efficiency the lattice map has the advantage of preserving moments, i.e. it does not introduce 
spurious diffusion into the problem. 

To test the features of the PDF seen in Fig. Q), we set 5 = mm(q s )/At. Note that as the advective map has 
no time scale, we set At = I to match an iteration of the mapping. Numerically this amounts to incrementing 
the mixing ratio at every site on the lattice by min(g s ) at every iteration. Now as min(g s ) is a function of a, 
6 decreases with increasing a. Interestingly, we observe that even though the maximum attainable mixing ratio 
is max((j s ), the system achieves equilibrium for much smaller maxima in q field. A snapshot of the equilibrium 
condensable field for a — 1,1.25 and 1.5 is shown in Fig. @ — note the increase in "graininess" of the field 
with a. The PDF's for these cases can be seen in Fig. (0 — the plotted PDF's are averages over the last 
couple of iterations of the map. In spite of the non-local nature of the mixing protocol the shapes of the PDF's 
display the features anticipated from the eddy-diffusive case, in particular we see the decrease in the slope of 
the PDF with S (the advective map is unchanged so K e is the same in all the simulations) also the character- 
istic roll-off induced via the form of the saturation profile in conjunction with the limit of rapid condensation is evident . 

Another measure of quantifying intermittency in a field is to examine its structure functions, defined as 
■Sn(H) =< \q(x + r) — q{x)\ n > (in the present situation < • > denotes a spatial average) higher moments (i.e. larger 
n) of the structure functions are most sensitive to the "roughest" regions of the field. For the condensable field, we 
observe that S n {r) ~ r^" for l\ < r < l 2 where l\ — > as we are dealing with a non-diffusive problem and l 2 is an 
outer scale which is smaller than the size of the domain. More to the point, the exponents Cn are anomalous in the 
sense Cn < ^Ci f° r n > 1- I n fact, we observe near saturation of the scaling exponents (i.e. £„ tends to a constant) 
for large n. Both the scaling of the structure functions and the extracted scaling exponents are shown in Fig. 
The anomalous behaviour is physically anticipated as smooth advection by itself leads to the formation of sharp 
fronts in finite time — hence one has a field composed of smooth regions interrupted by step like discontinuities. 
Indeed the combination of these two structures yields a profile which increases linearly for n <= 1 and then 
saturates to a constant (i.e. extreme intermittency) for n > 1 |24|. As the additional presence of condensation 
does not involve any smoothening of the field, we do not expect it to be able to mollify the discontinuities created 
via advection. But rapid condensation provides a bound for the magnitude of the jump across the advective 
discontinuity, i.e. \q(x + r) — q(x)\ < 1 — exp(— aL/2), hence we expect a weak dependence of Cn on a. Indeed, this 
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qualitative reasoning is bourne out in the lower panel Fig. In a similar vein, it has been shown (under further 
assumptions regarding the nature of the velocity field) that an interplay of advection with linear damping produces 
severe anomalous scaling Interestingly, such anomalous behaviour — with saturation of for large n — has 
been observed in an analysis of specific humidity fluctuations from high resolution aircraft data in the troposphere |'2o| . 



IV. CONCLUSIONS AND ATMOSPHERIC DATA 



We have examined the state of equilibrium achieved via the interplay of smooth advection and condensation. 
Exploiting the attainment of equilibrium without the presence of gradient fields allows us to make progress 
on the equation governing the PDF of the condensable substance. In particular, the limit of rapid condensa- 
tion admits a simplified Liouville equation with an integral normalization constraint on the PDF. Assuming a 
straightforward uniform solution to the governing equation shows the PDF to be tied to the particular saturation 
profile under consideration. Relaxing this uniformity and assuming an eddy-diffusive nature for the advective 
velocity fields shows the form of the PDF to still be controlled by the saturation profile but details, such as the 
rate of decay, to be sensetive to the eddy-diffusivity and strength of the forcing. Indeed, numerical simulations 
employing a lattice map for advective purposes reproduce these features. Further, scaling exponents extracted 
from structure functions of the condensable field show anomalous behaviour. Physically, the anomalous scaling is 
anticipated given the tendency of undiffused smooth advection to create sharp fronts. In fact, the near saturation 
of C„ for large n is a reflection of these fronts providing the leading contribution to the higher order structure functions. 

As mentioned in the Introduction, our motivation is to understand the distribution of water vapor along midlatitude 
isentropic surfaces. With this in mind we construct PDF's of the midlatitudinal specific humidity field along the 
300 K isentrope using data from the ECMWF reanalysis (ERA40) project. As is seen in Fig. (J5J, despite the 
simplifying assumptions made in our derivation, the PDF's from data can be approximated by the eddy-diffusive 
PDF. In fact, in Fig. (J5J we've also plotted the PDF that results from assuming uniformity (for the same a) and 
as is seen even though the primary shape of the data PDF follows from the saturation profile in conjunction with 
rapid condensation, it is the enhanced slope via ^f- that captures the decay of the PDF. Note that for large q the 
PDF's from data deviate significantly from the eddy-diffusive estimate. As we have considered the nondiffusive 
limit a possible source of this discrepancy might lie in the homogenization induced by diffusion. Secondly, as the 
source of water vapor lies in the tropics a boundary forcing might be more appropriate for the actual atmospheric 
problem. Even so, these moderately encouraging results lead us to conjecture that the ADC model — in the 
limit of rapid condensation with the proper saturation profile — driven by idealized velocity fields might be of 
use in predicting the statistical properties of the large scale distribution of water vapor in the midlatitude troposphere. 
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FIG. 1: Upper Panel : The PDF for the uniform case showing the dependence on a. Lower Panel : The PDF in i\'M) as a 
function of ^f- with a — 3. 




FIG. 2: A snapshot of the condensable field under the action of the map specified in l)14jl . The snapshots are after the PDF 
has settled into an invariant shape. The upper, middle and lower panels correspond to a = 1, 1.25 and 1.5 respectively with 
8 — min (q s ). 
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FIG. 3: The PDF's of the condensable field for a — 1, 1.25 and 1.5. The curves have been shifted for clarity. Note, as per Fig. 
Q, as a increases 5 decreases resulting in an increase in the slope (as n e is the same in every case) of the PDF. Also notice 
the increase in the sharpness of the peak with a. 
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FIG. 4: Upper Panel : Plots of \og(S q (r)) Vs. log(r) for the first six moments of the equilibrium condensable field with a = 2. 
Lower panel : Scaling exponents extracted from above (for a = 2) and similar plots (not shown) for a = 1, 1.5. 
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FIG. 5: The normalized specific humidity PDF's in the midlatitudes (between 30° and 60° in both hemispheres) along the 
300 K isentropic surface for Jan 97-01 from ECMWF data. The solid line with dots is a plot of 11311 (shifted vertically by a 
constant) with Sf- = 45, a = 3 and for comparison the dash-dash line is a plot of the PDF resulting from uniformity, again 
with a = 3. 



